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Abstract. A steady shear flow can drive supercooled liquids into a non-equilibrium state. Using molecular 
dynamics simulations under steady shear flow superimposed with oscillatory shear strain for a probe, 
non-equilibrium mechanical responses are studied for a model supercooled liquid composed of binary soft 
spheres. We found that even in the strongly sheared situation, the supercooled liquid exhibits surprisingly 
isotropic responses to oscillating shear strains applied in three different components of the strain tensor. 
Based on this isotropic feature, we successfully constructed a simple two-mode Maxwell model that can 
capture the key features of the storage and loss moduli, even for highly non-equilibrium state. Furthermore, 
we examined the correlation functions of the shear stress fluctuations, which also exhibit isotropic relaxation 
behaviors in the sheared non-equilibrium situation. In contrast to the isotropic features, the supercooled 
liquid additionally demonstrates anisotropies in both its responses and its correlations to the shear stress 
fluctuations. Using the constitutive equation (a two-mode Maxwell model), we demonstrated that the 
anisotropic responses are caused by the coupling between the oscillating strain and the driving shear flow. 
Due to these anisotropic responses and fluctuations, the violation of the fluctuation-dissipation theorem 
(FDT) is distinct for different components. We measured the magnitude of this violation in terms of the 
effective temperature. It was demonstrated that the effective temperature is notably different between 
different components, which indicates that a simple scalar mapping, such as the concept of an effective 
temperature, oversimplifies the true nature of supercooled liquids under shear flow. An understanding 
of the mechanism of isotropies and anisotropies in the responses and fluctuations will lead to a better 
appreciation of these violations of the FDT, as well as certain consequent modifications to the concept of 
an effective temperature. 

°: 

PACS. 05.70.Ln Non-equilibrium and irreversible thermodynamics - 61.43.Fs Glasses - 83. 50. Ax Steady 
shear flows, viscometric flow - 83.60.Df Nonlinear viscoelasticity 
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1 Introduction state. Even in strongly sheared non-equilibrium states, the 

structure and relaxation dynamics captured via the two- 
A comprehensive theory of systems driven into non-equilibriuBP 1 ^ correlation function exhibit very little anisotropy 
states is still under construction, in contrast to the well- CHI- This observation is in marked contrast to observa- 
established descriptions of equilibrium systems. Non-equilibritiffF of other complex fluids, such aspolymer solutions 
; states are generally characterized by violations of the fluctuate non-dense colloidal suspensions [10], m which struc- 
dissipation theorem (FDT). The FDT relates the response tural changes or anisotropic dynamics are induced by a 
functions to the associated correlation functions and holds drivin S shear flow ■ Furthermore, in glassy systems, 

in equilibrium states but is typically violated for non- including supercooled liquids, it has been suggested that 
equilibrium states. Much work has been devoted to under- th e equilibrium form of the FDT holds at long times with 
standing the relationship between response functions and the temperature T replaced by an effective temperature 
correlation functions in non-equilibrium situations; how- T ° ff OS' which indicates that T cS can be used to relate 
ever, this relationship remains unclear [THB]. the response and correlation functions. Several numerical 

It has been reported that supercooled liquids exhibit and theoretical works have examined the validity and the 
simple features even in non-equilibrium states. A steady role of the effective temperature in such situations UMI]. 
shear flow can drive supercooled liquids into a non-equilibrium Motivated by the above reports regarding the simple 

non-equilibrium properties of supercooled liquids, we in- 

a E-mail: e-mail: h-mizuno@cheme.kyoto-u.ac.jp vestigated the mechanical responses and the shear stress 

b E-mail: e-mail: ryoichi@cheme.kyoto-u.ac.jp fluctuations of a supercooled liquid in a non-equilibrium 
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state by means of molecular dynamics (MD) simulations. 
We first drove the supercooled liquid into a non-equilibrium 
state by applying a steady shear flow, and we then exam- 
ined the shear stress responses to oscillating shear strains 
in the sheared non-equilibrium state. In this study, we 
considered not only the weakly sheared situation but also 
the strongly sheared situation. In addition to the mechan- 
ical responses, the correlation functions of the shear stress 
fluctuations were also investigated in the non-equilibrium 
situation. We demonstrated the violation of the FDT and 
measured the magnitude of this violation using the effec- 
tive temperature. Shear stress responses and fluctuations 
are often useful for investigating non-equilibrium statisti- 
cal mechanics [6] . The aim of this study was to reveal the 
behaviors of the shear stress responses and fluctuations of 
the supercooled liquid in the non-equilibrium state. 

Several theoretical approaches addressing the mechan- 
ical responses of glassy systems are noteworthy, including 
the soft glassy rheology model [22J , the shear-transformation- 
zone theory [23], and the mode-coupling theory [23H25] . 
In Ref. [26] , the superposition rheology of glassy materials 
was investigated using the mode-coupling theory. Further- 
more, in the field of complex fluid rheology, several exper- 
imental studies have examined the mechanical properties 
of polymer solutions under a steady shear flow [27H3T] . 
More recently, Ref. |32] performed such an experimental 
study for glassy materials. We also note that constitu- 
tive equations, which detail the relationships between the 
stress tensor and the strain tensor, play an important role 
in predicting the fluid dynamics or the transport phenom- 
ena of the materials |33II34| . Several constitutive equations 
have been proposed to describe the mechanical properties 
of polymers 33. 3 1]. To the best of our knowledge, con- 
stitutive equations for supercooled liquids (or glassy sys- 
tems) have not yet been proposed for general shear strains 
(in tensor form) . In this study, we attempted to construct 
a constitutive equation for supercooled liquids that de- 
scribes our simulation results. 

The present paper is organized as follows. In Sect. [2] we 
briefly review our MD simulation. We also describe how to 
apply a steady shear flow and an oscillating shear strain. 
In Sect. [3J and 0] the results of the mechanical responses 
and the stress fluctuations are presented. In Sect. [3J we 
first indicate the mechanical responses obtained from the 
MD simulations. In this section, we also present a constitu- 
tive equation to describe our simulation results. In Sect. 2] 
we next show the results of the correlation functions of the 
shear stress fluctuations. Furthermore, we demonstrate 
the violation of the FDT and present the effective tem- 
perature as a metric for measuring the magnitude of this 
violation. Finally, in Sect. [5] we summarize our results. 



2 Simulation method 
2.1 Simulation model 

In this work, we performed MD simulations in three di- 
mensions. Our model liquid is a mixture of two atomic 



species, I and 2, with Ni — N2 — 5, 000 particles. The par- 
ticles interact via a soft-sphere potential <j)(r) = e(a a b/r) 12 
with (7 a i, = {a a + <7&)/2, where r is the distance between 
two particles, a a is the particle size, and a,ii G 1,2. The 
interaction was truncated at r = 3c a b- The mass ratio 
was set as m^/mi — 2, and the size ratio was set as 
02/01 — 1.2 to avoid system crystallization. Distances, 
times, and temperatures were measured in units of a%, 
t = (niiaf/e) 1 / 2 , and e/fcs, respectively. The particle 
density was fixed at a value of p — (Ni + Ni)/V = 0.8. 
The temperature was set to be T = 0.352 — 0.267. In this 
study, we mainly considered the state at the temperature 
T = 0.306. Note that the freezing point of the correspond- 
ing one-component model is approximately T = 0.772 |35) . 
At T = 0.352 — 0.267, the system is in a supercooled liq- 
uid state. After the system was carefully equilibrated un- 
der the canonical conditions, we applied a steady shear 
flow and an oscillating shear strain on the system using 
■the Lees-Edwards boundary condition [33] . We integrated 
the SLLOD equations of motion with the Lees-Edwards 
boundary condition, and the temperature was maintained 
by the Gaussian constraint thermostat 36 . The details 
of this simulation model can be found in previous stud- 
ies 0E] . 

2.2 Steady shear flow 

As mentioned above, after the quiescent equilibrium state 
was established, we applied a steady shear flow and an os- 
cillating shear strain. A steady shear flow is first applied to 
drive the supercooled model liquid into a non-equilibrium 
state [?H9] . We orient the x and y axes along the flow direc- 
tion and the velocity gradient direction of the steady shear 
flow, respectively, as shown in Fig.Q] We denote the shear 
rate of the steady shear flow as 7 SS , where the subscript 
"ss" indicates "Steady Shear how" . Figure [2] illustrates 
the shear rate 7 SS dependence of the shear viscosity 77 at 
various temperatures T = 0.352 — 0.267. The value of 77 
decreases with increasing 7 SS as rj ~ 7^" with v < 1.0, 
as demonstrated in previous studies [THS] . As observed in 
Fig. [21 the viscosity r\ displays a good fit to the functional 
form 77 = r} s o/(l + p^s) + Vf (the viscosity of constitutive 
Eq. ([6])). Note that this form 77 — r] s o/(l + fJ-jgs) + Vf i s 
the same as that proposed for pseudoplastic systems in 
Ref. [37]. 

In the present study, we mainly focused on two sheared 
non-equilibrium states, as indicated by the black circles in 
Fig. [2] One state occurs at T = 0.306 and 7 SS = 10 -4 , for 
which the shear flow is weak and thus the supercooled liq- 
uid is nearly in a Newtonian regime, i.e., the supercooled 
liquid is in the weakly sheared state. The other state con- 
sidered is at T = 0.306 and 7 SS = 10~ 2 , for which the shear 
flow is so strong that marked shear thinning occurs, i.e., 
the supercooled liquid is in the strongly sheared state. We 
note that a recent study [35] discussed the crossover from 
a Newtonian regime to a non-Newtonian regime (shear 
thinning regime) for sheared glassy systems in detail. To 
ensure that our simulations incorporated a different tem- 
perature case, we also considered the state at T = 0.267 
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10 3 , as indicated by the black square in Fig. 00]. We can calculate Gj? (uj) and G" y (uj) as the Fourier 

transformations of 8a lJ (t): 



2.3 Oscillating shear strain 

After the steady sheared state was achieved, we next added 
an oscillating shear strain to the main drive, as shown in 
Fig. [T] The oscillating shear strain was applied in a sinu- 
soidal form via the SLLOD algorithm [36]. We represent 
the oscillating shear strain as 8j t3 and its shear stress re- 
sponse as 8o %3 , where 7 1 - 7 and <j %3 are the ij components of 
the strain tensor 7 and the stress tensor <r, respectively, 
and ij — xy, xz, and yz. The difference in a quantity 
from its value in the absence of an oscillating strain is de- 
noted as 8. Notably, even in the absence of an oscillating 
strain, a xy has a value a ss = ?7(7ss)7ss due to the steady 
shear flow 7 SS , and we should therefore calculate 8a xy as 

foxy = a xy _ ^ ( 6a xz = Q xz &nd S(J yz = a yz _) jjg^ 

we stress that according to the three components ij = xy, 
xz, and yz, there are three different ways to apply an os- 
cillating shear strain, as shown in Fig. [T] Figs. [Ha), (b), 
and (c) correspond to ij — xy, xz, and yz, respectively. 
The strain tensor, 7 = Vit + *(Vw), is written as Eqs. ([T]), 
©, and ([3]) for ij — xy, xz, and yz, respectively. 
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7 = 



7 



KstSY* 

j ss t 
5~f xz 

j ss t 

i ss t w 
8-/ yz 



(1) 



(2) 



(3) 



In the present situation, in which a steady shear flow with 
7 S;5 is applied, the three stress responses, 8a xy to 8j xy , 
S(T XZ to 5~f xz , and 8a yz to 8j yz , are generally different, 
whereas these three responses are exactly the same in the 
equilibrium state 7 SS = 0. 

The oscillating shear strain was expressed in a sinu- 
soidal form: 

8^3 =7 sin(W), (4) 

where 70 and uj are the amplitude and the frequency of 
the oscillating strain, respectively. In this study, we set the 
amplitude 70 to be 70 = 0.01 — 0.2. The amplitude 70 = 
0.01 is small enough that the response 8a %3 is linear with 
respect to Sj 10 . In contrast, if 70 > 0.01, 8<j %3 becomes 
non-linear with respect to 8j Z3 . 

It is beneficial to use the shear moduli G^(uj) and 

G^(uj) instead of the full time history 8a 13 (t), where 
the subscript "7 SS " denotes a value in the sheared non- 
equilibrium state. The values G 1 ? 3 (uj) and G 1 ^ 3 (uj) are the 
storage modulus for the elasticity and the loss modulus 
for the viscosity, respectively, and are often used to mea- 
sure the viscoelastic properties of the materials 24,25,39, 



uj f n ' u 8a 13 (t) 



7r/u 



7o 



u r /u 8a 13 (t) 



7r/u 



7o 



sm(ujt)dt, 



cos(ujt)dt. 



(5) 



In this case, the time history 8a %3 (t) can be expressed 
as 8a i3 (t) = G^f a (o;)n)Bin(wi) + G'^u)^ oos(wt). The 

shear moduli G^(ui) and G"^(w) depend on the three 
quantities uj, 70, and 7 SS . If the amplitude 70 of the oscil- 
lating strain is small and the steady shear flow 7 SS is weak, 
then G'? 3 (uj) and G" u (uj) reduce to the linear shear mod- 
uli and depend only on the frequency uj of the oscillating 
strain. However, when the amplitude 70 becomes large or 
the steady shear flow 7 SS becomes strong, then signifi- 
cant non-linearity arises, so Gl^ (uj) and G"^ (uj) become 
the non-linear shear moduli [21,25,39,40], which depend 
not only on the frequency uj but also the amplitude 70 
or the steady shear rate j ss . We should note that when 
the responses are non-linear with respect to the oscillat- 
ing strain, there are higher harmonic contributions to the 
responses [25][3j3][30] . In the present study, we considered 
only the first harmonic contribution, i.e., the shear mod- 
uli G? 3 s (w) and G"^(w), which were verified to exert a 
dominant effect compared with the higher harmonic con- 
tributions. 



3 Result I: Mechanical responses 

In this section, the results of mechanical responses are 
discussed. We present the shear moduli and G^ 

obtained from MD simulations. The simulation cases are 
summarized in Table [1] As we mentioned in Sect. [2~2l we 
primarily focused on two sheared states with the temper- 
ature T = 0.306: the weakly sheared state 7^ = 10 -4 and 
the strongly sheared state 7 SS = 10~ 2 . In addition to these 
two states, we considered another state with T — 0.267 as 
a case involving a different temperature. We also present 
the constitutive equation that captures the key features of 
the simulation results. 



3.1 Mechanical responses in the weakly sheared state 

We first present the results of the mechanical responses 
in the weakly sheared state T = 0.306 and 7^ = 10~ 4 . 
Figure [3] illustrates the shear moduli G'V and G"^ at the 
small amplitude 70 = 0.01 of the oscillating shear strain. 
The value 70 = 0.01 is small enough that the mechanical 
response is linear with respect to the oscillating strain. 
In the same figure, we also show the values of the shear 
moduli G' and G" in the equilibrium state to clarify the 



effects caused by the steady shear flow 7 SS . As can be seen 
from Fig. [3J the shear moduli G' eq and G" q demonstrate 
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the typical dependence on frequency uj of the Maxwell 
model with two time scales [33l[34] . In Fig. [3l we indicate 
these two time scales as t s q (slower time) and 77 (faster 
time). As is well known, the stress correlation functions of 
supercooled liquids exhibit two-step relaxation [4TH2"] (see 
also Fig. IH]). The slower relaxation is called cv-relaxation, 
and its relaxation time is thus known as the cv-relaxation 
time r Q . The faster relaxation results from the thermal 
vibrations of the particles, and its time scale is known as 
the Einstein period te (uJe = 27T /te is the Einstein fre- 
quency). The two time scales r s o and 77 are equivalent to 
r a and te, respectively (t s q ~ r Q 3> 77 ~ te)- Here, we 
note that at approximately the slower time scale t s q ~ r a , 
supercooled liquids exhibit the crossover from liquid-like 
behavior to solid-like behavior. At low frequencies uj, G" q 
is larger than G' eq , which results in liquid-like behavior. 
As uj becomes large, G' eq becomes larger than G" , which 

results in solid-like behavior. At uj ~ , G' cq ~ G" q , i.e., 
the crossover occurs. In addition, as we will explain in de- 
tail in Sect. 13.31 and G-^ can also be well described 
by the two-mode Maxwell model as in Eq. (O; therefore, 
the mechanical properties of the supercooled liquid can be 
characterized by two time scales not only in the equilib- 
rium state but also in the sheared non-equilibrium state. 
In Fig. [31 two time scales in the sheared state are indicated 
as t s (slower time) and 77 (faster time). 



As in Fig. [21 although we can recognize the effects 
due to the steady shear flow at the very low frequencies 
of uj < t~ , these effects are very small. In the whole 
frequency uj region except for oj < r s _1 , all and G"^ 
values almost coincide with G' eq and G" q . The slower time 
scale t s (equivalent to r a ) is also close to the equilibrium 
time scale of t s q, although t s is a little shorter than t s q. 
The faster time scale (equivalent to te) is unchanged at 
77. Thus, the weak steady shear flow with 7 SS = 10~ 4 
produces only small effects on the mechanical responses 
in the very low frequency region uj < t~ 1 . 



Furthermore, in Fig.[H we demonstrate and G" y 
at the large amplitude 70 = 0.1 of the oscillating shear 
strain. The value 70 = 0.1 is large; thus, the mechanical re- 
sponse is non-linear with respect to the oscillating strain. 
Comparing Fig. H ( 7o = 0.1) with Fig. [3j ( 7o = 0.01), we 
can observe that the storage modulus decreases for a 

greater 70, and the loss modulus G" lJ becomes larger rela- 
tively. This observation implies that the larger oscillating 
strain makes the supercooled liquid more liquid-like. Such 
non-linear viscoelasticity has also been observed in soft 
materials [2~4"ll3l?] , dense colloidal suspensions 25 , and su- 
percooled polymer melts [10] . In addition, it is also notable 
that the effects of the steady shear flow j ss become smaller 
for a greater 70 = 0.1, as the oscillating shear strain be- 
comes relatively strong compared with the steady shear 
flow. 



3.2 Mechanical responses in the strongly sheared state 

Figures [5j and [HI illustrate the results for GJ? and G"^ 3 in 

the strongly sheared state T — 0.306 and j ss = 10~ 2 . The 
amplitude 70 of the oscillating shear strain is 70 = 0.01 
(linear regime) in Fig. [3] and 70 = 0.1 (non-linear regime) 
in Fig. [5] As in Figs. [5] and [Bl we can easily recognize 
the effects resulting from the steady shear flow at low fre- 
quencies uj < t^ 1 for all and G"*^, whereas at high 
frequencies uj > tJ 1 , these effects are not observed, i.e., 

all G'? J and G" y values coincide with G' cq and G" q . No- 
tice that due to the steady shear flow, the slower time 
scale (equivalent to t q ) becomes dramatically shortened 
from t s q to t s , whereas the faster time scale (equivalent 
to te) is unchanged at Tf. The strong shear flow with 
7 SS = 10~ 2 influences the mechanical responses to a much 
greater extent than is observed for the weak shear flow 
with 7 SS = 10~ 4 . 

From Figs. [5Jand[Hl we obtain two remarkable results 
for both the linear and the non-linear responses (both 70 = 
0.01 and 70 = 0.1). First, two components, xz and yz, of 
G? 3 and G" y coincide with each other surprisingly well, 
even at low frequencies oj < r" 1 (refer to the upper and 
lower triangles in Figs. [5] and [6]). Despite the strong steady 
shear flow, two of the stress responses, 8a xz to 8^f xz and 
Stjyz to S"/ yz , are the same. This result demonstrates the 
isotropic aspect of this system. Second, the behavior of the 
xy component (refer to the circles in Figs. [5] and [5]) is quite 
different from those of the xz and yz components. The xy 
component is smaller than either the xz or yz component 
at low uj. In addition, as uj decreases, the storage modulus 
G^ decreases much more rapidly than G^ and G . At 
low uj, G^ takes on negative values (data are not shown 
in Figs. [S] and [5]). Thus, due to the steady shear flow, the 
mechanical properties of the xy component are notably 
different from those of the other components. This result 
demonstrates the anisotropic responses of this system. We 
will discuss the origin of the anisotropic responses, i.e., the 
difference between the xy component and the xz and yz 
components, in Sect. 13.41 Similar behaviors of G'* v and 
G"^ v have been previously observed in polymer solutions 

ESEo]. 

In Figs. [7] and we present and G"^ at the 

different temperature case of T = 0.267 and j ss = 10 -3 , 
for which the supercooled liquid is also in the strongly 
sheared state. Figures [7] and [5] show the results for the 
different amplitudes 70 = 0.05 and 0.2, respectively. We 
can confirm that the same observations from Figs. [5] and [6] 
are also obtained in the scenario displayed in Figs. [7] and 

El 



3.3 Constitutive equation 

We attempted to construct a constitutive equation de- 
scribing our simulation results and obtained the following 
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two-mode Maxwell model equation: 



da %3 

r\ +r s (7)-jr =v s {i)i lJ , 



dt 
da % l 

°t + T f-£- = vfy v > 



(6) 



with 



r) s (i) = 



TsO 



1 + fi(j x y 2 + j xz2 + 'jy z2 ) u / 2 ' 

VsO 

i + n^y 2 + j xz2 + ^v z2 y/ 2 ' 



(7) 



where ij — xy, xz, or yz. The model equation consists 
of a slower component and a faster component, which 
are denoted by the subscripts "s" and "/", respectively. 
The stress tensor cr 4 - 7 is written as cr u = a l J + aj . As 
in Figs. [3] or [5] the (linear) shear moduli G' eq and G" q in 
the equilibrium state demonstrate the typical frequency 
ui dependence of the Maxwell model with two character- 
istic times [331134"] ; therefore, we considered the two- mode 
Maxwell model equation. In this model, the shear viscosity 
rj and the instantaneous shear modulus Goo are described 
as 7] = r) s +r)f and Goo = G s + G f = Vs/T s +Vf/ T f> respec- 
tively. The values rj s and G s — t] s /t s are slower compo- 
nents of 77 and Goo, whereas the values 77/ and Gf = 77/ /r/ 
are faster components. As we mentioned in Sect. 13.11 the 
two time scales t s and 77 are considered to be equivalent 
to the a-relaxation time T a and the Einstein period te, 
respectively. Therefore, we naturally assumed that the val- 
ues Tf and rjf characterizing the faster component are con- 
stant and unaffected by the applied shear flows, whereas 
the slower components r s and rj s do depend on the ap- 
plied shear flows. We set t s and rj s as functions of only 
the total strength of the shear rate (j xy2 + i xy2 + i' xy2 ) 1/2 
as in Eq. ([7]), which reflects the isotropic feature observed 
in MD simulations, i.e., that the shear responses of the 
xz and yz components coincide with each other in Figs. 
[5] and [6] In addition, we assumed r s ~ rj s and set the 
same functional form for t s and t] s as in Eq. (J7J. The 
functional form simply arises from the shear-thinning be- 
havior shown in Fig. [2] In fact, when we consider a steady 
shear flow with 7 SS , the shear viscosity 77 is described as 
V = r) s + r\f = 77 s0 /(l + yrfss) + Vf, which is precisely 
the shear-thinning form. As we mentioned previously, this 
functional form of 77 is the same as that proposed for pseu- 
doplastic systems in Ref. [37] , 

Together, the constitutive Eq. © and Eq. (J7J have six 
parameters, 77^, ??/, G s0 = Vso/t s o, G f = 77//V/, /j, and 
v. Four of these parameters, 77 s o, ??/, G s o, and G/, charac- 
terize the linear mechanical responses in the equilibrium 
state. We can determine these four parameters from the 
shear stress correlation function G(t), defined as 



(8) 



cq 



where 5a represents the shear stress fluctuations and () 
denotes the ensemble average in the equilibrium state. The 
shear viscosity rj and the instantaneous shear modulus G c 



are related to the function G(t) through 77 = J °° G(t)dt 
and Goo = G(t = 0) [43]. In addition, as is well demon- 
strated in Ref. [U] and Fig. |H1 the slower relaxation of 
G(t) can be well fitted by the stretch exponential form 
G Q exp (—(t/T a )^), where the value G a is known as the 
plateau modulus . We therefore determined 77^0 and 

G s o (the slower components of the shear viscosity and 
modulus) as 



VsO 



G a exp(-(t/r a f)dt : 



(9) 



G s o — G a , 



and the values of 77/ and Gf (the faster components) were 
determined as 

i* 00 

77/ = 77 - -qso = / G{t) - G a exp (-(i/r Q ) v ') dt, 

Jo UUJ 

Gf = Goo — G s o = G(t = 0) — G a . 

Figure ITT1 shows the temperature T dependences of 77^0 , 
77/, G s o, and Gf. The slower component 77,50 increases 
dramatically with decreasing temperature T The 
faster component 77/ is much (several orders of magnitude) 
smaller than ?7 s o, and the shear viscosity 77 is almost the 
same as the slower component 77 s o. On the other hand, the 
values G s o and G/ are nearly constant with respect to T, 
as was previously observed in Ref. [4"4"] . 

The remaining two parameters, fi and v, characterize 
the non-linearity resulting from the driving shear flow and 
can be determined from fitting the shear viscosity 77 = 
??so/(l+At7ss)+ 7 ?/ of the model equation to the simulation 
data ?7(7 SS ) shown in Fig. [5J The functional form 77 = 
Tjso/0- + M7ss) + Vf i s transformed to 



1 



77-77/ 



1 

VsO 



JL 

VsO 



(11) 



and thus, we fitted the function {^/Vso)iss to the data 
l/(r)(iss) — Vf) ~ l/^so, as in Ref. j37j. As shown in Fig. 
HU1 the straight line is well fitted in the log-log plot, for 
which we performed the least-squares fit. The slope and 
intercept of the fitted line then correspond to the values 
v and n/vsQ, respectively. As can be observed in Fig. [5J 
the function 77 = 7? s o/(l + MTss) + Vf with the obtained 
values of /i and v demonstrate a good fit to the simula- 
tion data 77(7^5). FigurefTTIalso illustrates the temperature 
T dependences of /i and v. The value /i increases drasti- 
cally with decreasing T in a similar way as the viscosity 
77, whereas v is insensitive to T and takes values between 
0.8 and 1.0 [7J{5]. 

We calculated the shear moduli and G"*^ using 
the constitutive equation ^ with the parameters pre- 
sented in Fig. [IT] Figures [3] and [4] also show the results of 
G^ s and G obtained from Eq. (O in the weakly sheared 
state of T — 0.306 and ^ ss = 10~ 4 . By comparing the 
lines (constitutive equation) with the symbols (MD sim- 
ulation), we can clearly observe that Eq. ([6]) reflects the 
results of the MD simulations quite well in both the equi- 
librium state and the weakly sheared state. Notice that the 
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constitutive equation can account for not only the linear with 
responses (the amplitude 70 = 0.01 in Fig. [3]) but also the 
non-linear responses (70 = 0.1 in Fig. 2]). In addition, we 
show G]? and G"^ 3 of Eq. (0 in the strongly sheared state 
T = 0.306 and 7,**= 10~ 2 in Figs. [5] and © These results 
demonstrate that the constitutive equation also functions 
surprisingly well even in the strongly sheared state, ex- 
cept for the storage modulus Gj* v at low frequencies oj. 



The modulus G'? v takes on negative values at low oj, as 
we mentioned previously in Sect. 13.21 and Eq. © cannot 
account for this negative storage modulus. At this stage, 
we do not understand the mechanism and the origin of the 
negative values in question for Gj^ , and this topic will be 
a subject of future work. Once the mechanisms underlying 
these values are more fully elucidated, we will be able to 
propose certain modifications of the constitutive equation 
to account for this negative modulus. 

Furthermore, we verified the validity of the constitu- 
tive equation in a case involving a different temperature 
T = 0.267 and 7 SS = 10~ 3 , as shown in Figs. and E 
In Figs. [7] and [51 we once again observe that the consti- 
tutive equation is valid except for the storage modulus 
G'? v at low frequencies oj. We stress that all six parame- 
ters of the constitutive equation have physical significance 
and can only be completely determined by the mechanical 
properties (G(t) and r]('j ss )) in the equilibrium situation 
and the steady sheared situation. Using these six parame- 
ters, we can accurately predict mechanical properties in 
more general situations, e.g., under two different shear 
strains (the steady shear flow and the oscillating strain). 
The constitutive Eq. (jB]) is much simpler than other equa- 
tions obtained for typical complex fluids such as polymer 
solutions |33ll34j . and interestingly, even in the strongly 
sheared state, the mechanical responses can be well fit by 
this simple constitutive equation. 



3.4 Origin of anisotropic mechanical responses 

As demonstrated in Figs. [5] and [6] (and also in Figs. [7] and 
[FJ) , G'^ s and G"^ differ from the other xz and yz compo- 
nents at low frequencies oj in the strongly sheared state. 
We can understand the origin of this difference through 
the examination of the constitutive equation ©. At low oj, 
the slower component a\ 3 is dominant compared with the 
faster component aj ; thus, in the following analysis, we 

consider only a l J . To reveal the effects due to the steady 
shear flow, we analyzed the constitutive equation under 
the condition 70^ <C j ss . This condition, 70a; <C j 38 , indi- 
cates that the shear rate of the steady shear flow is much 
larger than that of the oscillating strain, which holds true 
for low oj. Under this condition, the following two equa- 
tions were obtained: Eq. (|12p for the xy component and 
Eq. (|13p for the xz and yz components. 



dSa xy 



at 



W x \ (12) 



1 + M7ss 



Vs 



1 



(14) 



Here, 8j 13 — joojcosojt is the shear rate of the oscillating 
shear strain. Note that the equations for the xz and yz 
components are exactly identical, which arises from con- 
stitutive Eq. © being isotropic with respect to the strain 
tensor 7. By comparing Eqs. (jT2"j) and (fTS"]) . we see that 
the difference between the xy component and the other xz 
and yz components originates from the second term of the 
right-hand side of Eq. (JT2J, i.e., the "7^7^" term. This 
term arises from the coupling between the driving shear 
flow 7 SS and the oscillating strain E^ xy . In the xz and yz 
components, such a coupling does not appear because the 
driving shear flow and the oscillating strain impact sepa- 
rate components. However, in the xy component, the driv- 
ing shear flow ^ ss and the oscillating strain 5 / y xy affect the 
same components, and the coupling term :l jssdi xv " there- 
fore arises, which is the origin of the difference between the 
xy component and the other xz and yz components. 

To conclude this section, we allude to certain studies 
investigating the mechanical responses of polymer solu- 
tions under a steady shear flow [6l l27ti3T] . in which two 
responses 5a xy to 6j xy and 5a yz to Sj yz have been ex- 
amined. Note that in these studies, the responses Sa xy to 
8^f xy and 5a yz to 5^j yz are called "parallel superposition" 
and "orthogonal superposition" , respectively. It was re- 
ported that the steady shear flow causes an effect called 
the convective constraint release effect [3D] or induces an 
anisotropic mobility [5J, which both accelerate the time 
scale of polymer dynamics. This result is similar to our 
findings that the driving shear flow makes the character- 
istic times faster. However, in polymer solutions, the ori- 
gin of the difference between the xy and yz components 
is still not clear. Furthermore, to the best of our knowl- 
edge, there are no studies investigating the mechanical 
responses of polymer solutions for the xz component. It is 
expected that polymer solutions exhibit different mechan- 
ical responses in the xz and yz components, in contrast 
to supercooled liquids, which surprisingly exhibit the same 
responses for these components. 



4 Result II: Stress fluctuations 

In this section, the results of stress fluctuations are pre- 
sented. We show the stress correlation functions in the 
sheared non-equilibrium state. In addition, we also demon- 
strate the violation of the FDT and present the frequency- 
dependent effective temperature as a metric indicating the 
magnitude of this violation. 



4.1 Stress correlation function 

We examined the correlation function of the shear stress 
fluctuations, defined as 



(ij = xz, yz), (13) 



(t) = ^(Sa^(t)6a i3 (0)). g 



(15) 
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where 8a lJ represents the shear stress fluctuations in the 
ij component of the stress tensor <x and ()^ denotes the 
ensemble average in the sheared non-equilibrium state. 
The functions GJP, and are plotted in Fig. 

[T21 which also shows G cq in the equilibrium state. The 
temperature is T — 0.306. The shear rate 7 SS is 7 SS = 
10 -4 (weakly sheared state) in Fig. [T^Ta) and j ss = 10~ 2 
(strongly sheared state) in Fig.Q^b). In the weakly sheared 
state (Fig. rnZT a)). although G?* behaves slightly different 

from G eq , there are only small differences between G? 



and G cq ; thus, the weak steady shear flow with j s 



10" 



produces only small effects on the shear stress fluctua- 
tions, as was observed in the mechanical responses shown 
in Figs. |3] and [4] However, in the strongly sheared state 
(Fig. rnZT b)). we can easily recognize that all Gt? behave 
much differently than G eq . Due to the strong shear flow, 



cq- 



not yet understand the mechanism of either the negative 
fluctuation correlation or the negative shear modulus at 
this stage. 



4.2 Violation of fluctuation-dissipation theorem 

The relationship between the response and correlation func- 
tions is of great interest and importance [1T J6TT3TI2T] . If we 
use the shear moduli G? 3 and G" u as the response func- 
tions, the associated correlation functions are the Fourier 
transforms of the correlation function G\ 3 , defined as 



G 



(w) = u / CtP (t) sm{ut)dt, 



G^ r ,H=W Gt(t)cosM)cft, 



(16) 



all of the G l ? a values relax rapidly compared with G, 
It should be noted that even for a small time t < r f , all where the subscript w denotes the « corre l a tion func- 

° l L values differ from which means that thc stcad y tion". We note that the response functions must be G* 



shear flow with 7 SS = 10 -2 affects the fluctuations even 
at t < Tf ~ 0.15, time scales that are much smaller than 
--y^ 1 = 10 2 . A similar result has been observed in sheared 
foam [TBI. This trend is different than those of the re- 
sponses G? 3 and G"^ 3 shown in Figs. [5] and H2 for which 
at short time scales (u) > t7 ), all G'* 3 and G" y values 
coincide with G' eq and G" q , i.e., effects due to the steady 
shear flow are not observed. 

Fig- HHb) demonstrates that G?^ and G^ z are identi- 
cal, whereas G* v behaves quite differently and even takes 
on negative values. The coincidence of G?^ and GV^ demon- 
strates that the shear stress in the xz and yz components 
fluctuates in the same manner despite the strong shear 
flow in the xy component. Such isotropic fluctuations have 
also been observed in the density correlation function and 
the mean square displacement [7]-[9]. In contrast, the dif- 
ferent behavior of G^ y demonstrates an anisotropy, which 
has also been detected in the four-point correlation func- 
tion [47]. Similar isotropy and anisotropy were also ob- 
served in the mechanical responses and G"^ depicted 
in Figs. [5] and HJ As discussed for Eqs. (TT2J) and (TT3")) . the 
coupling between the driving shear flow j ss and the oscil- 
lating strain 8^ xv makes the response in the xy component 
different from those in the xz and yz components. We can 
expect that similar to the effect seen for the responses, 
coupling interactions between the driving shear flow and 
the fluctuations of the xy component also arise, whereas 
such coupling does not appear in the xz and yz compo- 
nents. 

In addition, it should be noted that the stress corre- 
lation function G^ does not decay monotonically and 
takes on negative values at long time scales in Fig. H2T b). 
The same relaxation behavior has been also reported in 
sheared foam fT^] . As we observed in Fig. [5] and HI the 
storage modulus G'* v also becomes negative at long time 
scales (low frequencies ui). We speculate that the negative 
shear modulus and the negative fluctuation correlation of 
the xy component may be interrelated. However, we do 



and G"^ 3 at the small amplitude 70 = 0.01 of the os- 
cillating strain, i.e., the linear responses. Figures 1131 and 
H4l present the correlation functions G c ^ r ^ and G c £ ■ 

with the response functions G'^ 3 and G"^ 3 in the weakly 
sheared state (T — 0.306 and ^ ss — 10 -4 ) and the strongly 
sheared state (T = 0.306 and j ss — 10 -2 ), respectively. 
The functions in the equilibrium state are shown in the 
same figures. We can observe that the response and cor- 
relation functions coincide with each other in the equi- 
librium state, which implies that the FDT holds. In con- 
trast, the response and correlation functions do not coin- 
cide in the non-equilibrium state, and therefore, the FDT 
is violated. The large violation of FDT is observed in the 
strongly sheared case (Fig. IT4"|) . Note that in the strongly 
sheared case, the FDT is violated even at high frequencies 
u> > tJ 1 for the storage modulus G'? 3 s and G^ r ^ ss . This 
fact can be more clearly observed by noting that the effec- 
tive temperature T e ' ff , measured from the storage modulus 
in Eq. (flTl) . has slightly higher values than the tempera- 
ture T = 0.306, even at high u > tJ 1 , as shown in Fig. 
[TBI The violation of the FDT is physically interpreted as 
entropy production [3d], a steady-state probability cur- 
rent [3], or other related physical quantities. Recent works 
have attempted to provide a unifying framework for this 
phenomenon [H[S]. 

4.3 Frequency-dependent effective temperature 

To examine the concept of an effective temperature [13i - 
I2T] , we defined the frequency-dependent effective temper- 
atures T' eS {uj) and T" s (uj) as 



») = Or7» 

T G ">) 



(17) 
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Note that the effective temperatures T' eS and T" s are de- 
fined from different observables, i.e., the storage modulus 
and the loss modulus, respectively. Figures [T5] and [H)] il- 
lustrate the frequency uj dependencies of T c ' ff and T" s in 
the weakly sheared state (T = 0.306 and 7 SS = 10 -4 ) and 
the strongly sheared state (T = 0.306 and 7 SS = 10~ 2 ), 
respectively. In the equilibrium state, T^ ff and T" s are ex- 
actly the same as the temperature T = 0.306 for all fre- 
quencies uj, as is well demonstrated in Figs. [15] and [TCI 
Thus, the temperature T exactly relates the response func- 
tion to the correlation function in the equilibrium state. 
However, in the non-equilibrium state, both T c ' ff and T" s 
differ from T = 0.306, particularly in the strongly sheared 
case (Fig. \W\i . If the concept of an effective temperature 
is valid, then, at low frequencies to (long time scales), T' eS 
and T" s should coincide with each other and have the same 
value for all components ij = xy, xz, and yz, i.e., one 
scalar quantity should relate the response function and the 
correlation function of any observable and any component. 
However, as shown in Figs. [TBI and [11)1 T ' ff and T" s differ 
from each other and have different values between compo- 
nents. Note that in the strongly sheared state (Fig. [TB"]) , 
T^g of the xy component has negative values at low fre- 
quencies uj because of the negative storage modulus G'? J 
(data are not shown in Fig. ITB"]) . The differences in the 
effective temperatures have been also observed between 
different observables [T5l[T7] and different directions [2"T] . 
This result indicates that the use of only one effective tem- 
perature (one scalar quantity) cannot completely charac- 
terize the relationship between the response functions and 
the correlation functions. For the supercooled liquid, it is 
notable that the effective temperatures T e ' ff and T" s have 
the same values for the xz and yz components over all 
frequencies uj. This isotropic feature is considered to be a 
characteristic of supercooled liquids (or glassy systems) . 

In addition, it is interesting that the frequency uj de- 
pendencies of the effective temperatures T' cS and T" s no- 
tably differ from each other over the entire domain of uj. 
The value T e 'g coincides with T = 0.306 at high frequen- 
cies uj > t ~ . As uj decreases, T" s starts to deviate from 
T = 0.306 and approaches higher values: 0.350 for the 
xz and yz components and 0.550 for the xy component 
in the weakly sheared case (Fig. [T5J and 0.450 for the 
xz and yz components and 0.800 for the xy component 
in the strongly sheared case (Fig. [TB"]) . This behavior of 
T" s is consistent with previous results [LH[T4l,[20l[^T] . i.e., 
the effective temperature is equivalent to the bath tem- 
perature at short times (high uj) and a higher tempera- 
ture at long times (low uj). However, the value T' cS can 
be slightly higher than T and does not coincide with T 
even at high uj > tJ 1 , as can be clearly observed in the 
strongly sheared case (Fig. [16]). As uj decreases, T e ' ff be- 
gins to increase at a frequency that is much faster than 
t" 1 (near uj — tJ 1 ~ t^ 1 in the strongly sheared case). In 
addition, Tg ff for the xy component does not approach a 
positive value. Even in the weakly sheared case (Fig. [15]) . 
we did not obtain the asymptotic value within our simula- 
tions. These results for T' eS are inconsistent with previous 
results 13,14.20,21]. Therefore, the effective temperature 



measured in previous studies [T3llT4l[2lIl[2T] may be physi- 
cally similar to T" s but meaningfully different from T' eS . 

To demonstrate that the effective temperature mea- 
sured in previous works [IS] [TH (SHUT] is quantitatively 
identical to the present T" s , we investigated the mechani- 
cal response to a small constant shear strain. At time t = 
0, the constant small shear strain Sj^ — 0.01 was applied 
to the sheared supercooled liquid. Note that the response 
is linear with respect to the shear strain Sj^ = 0.01. We 
calculated the response Sa^ using MD simulation and ob- 
tained the susceptibility x^ sa (*) as 



(S^(t)-6gv(0)). s 



(18) 



where () • denotes the ensemble average in the sheared 
non-equilibrium state. The effective temperature T e g is 
defined as the inverse slope of the plot of the susceptibility- 
correlation function for long times [13l[T4l[20l[21] , as given 
by Eq. (CET 



4L(*) 



1 



off 



(TG^ (O)-TGJ? (i)). 



(19) 



Figure [T7] shows the plot of the susceptibility-correlation 
function for the ij = xy, xz, and yz components. The 
shear rate 7 SS is j ss = 10~ 4 (weakly sheared state) in 
Fig. [TWa) and j ss = 10~ 2 (strongly sheared state) in Fig. 
[TTTb'). In the same figure, the plot of the susceptibility- 
correlation function in the equilibrium state is also shown. 
In the equilibrium state, the inverse slope is exactly T = 
0.306 in the whole region. In contrast, in the sheared non- 
equilibrium state, the inverse slopes for long times take 
values higher than T = 0.306, and these values coincide 
well with T" s {uj) at low frequencies uj: 0.550 for the xy 
component in the weakly sheared case (Fig. H7Ta)) and 
0.450 for the xz and yz components and 0.800 for the xy 
component in the strongly sheared case (Fig.H7Tb)). Note 
that in the weakly sheared case, the inverse slopes for the 
xz and yz components should be 0.350 (T" s (uj) at low 
uj), although we cannot distinguish such slopes numeri- 
cally. Thus, we can conclude that the effective tempera- 
ture T e g used in previous works P^l[T4l[^rjll^T] and T e 'g are 
exactly same. In fact, the effective temperatures T fj and 
T" s are mathematically connected, as the responses to a 
small constant strain and an oscillating strain are related 
by a Fourier transformation. 



5 Conclusion 

In this study, we examined the shear stress responses and 
fluctuations of a supercooled liquid in sheared non-equilibrium 
states. Interestingly, for the two components differing from 
that of the driving shear flow, the same responses and fluc- 
tuation correlations were observed despite the presence of 
that driving shear flow. From these responses and fluctu- 
ations, identical magnitudes of the violations of the FDT 
were obtained. These results demonstrate the isotropic as- 
pect of this system [7HS]- Based on this isotropic feature, 
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we successfully constructed the two mode-Maxwell model 
with isotropic non-linearity, written as Eq. ([6]) and Eq. 
(J7J|, for the constitutive equation of supercooled liquids. 
This simple constitutive equation is surprisingly accurate 
at describing the mechanical properties of the supercooled 
liquid not only in the quiescent equilibrium state but also 
in the sheared non-equilibrium state. 

In contrast, for the same component as that of the driv- 
ing shear flow, both the responses and fluctuations are 
quite different from those of the other two components. 
This result demonstrates the anisotropic aspect of the 
system [37]. Using the constitutive equation, we demon- 
strated that the anisotropy in the responses derives from 
the coupling between the steady shear flow and the oscil- 
lating shear strain, a coupling that only arises when both 
factors impact the same component. Coupling with the 
driving shear flow can also be expected for the stress fluc- 
tuations, which may cause anisotropic fluctuations, just 
as such a coupling caused anisotropic responses. This cou- 
pling between the driving shear flow and the stress fluctu- 
ations will be a future topic. Furthermore, for this compo- 
nent, we observed a negative storage modulus and a neg- 
ative fluctuation correlation, which may be interrelated 
phenomena. The origin and relationship between these 
negative values of the storage modulus and fluctuation 
correlation will also be addressed in future study. Due 
to its anisotropic responses and fluctuations, the magni- 
tude of this component's violation of the FDT differs from 
those observed for the other components [21] , which indi- 
cates that a simple scalar mapping, such as the concept 
of an effective temperature [T3H2T] , oversimplifies the true 
nature of supercooled liquids under shear flow. In fact, 
we quantified the magnitude of this violation using the 
frequency-dependent effective temperature, defined as Eq. 
(fTT)) . and demonstrated that the effective temperature ex- 
hibits notably different values for different components. 

Finally, we stress that even in the strongly sheared 
situation, supercooled liquids exhibit highly isotropic na- 
tures (isotropic dynamics and structure [7_HS] as well as 
isotropic mechanical responses and fluctuations), which 
cannot be expected for typical complex fluids in which 
the driving shear flow induces anisotropic dynamics or 
a structural change [ITJ[l2]. The constitutive equation, 
based on the isotropic features of supercooled liquids, is 
much simpler than those proposed for polymer solutions 
[27-31,33,34 . These isotropic and simple features are con- 
sidered to be characteristic of the supercooled liquids (or 
glassy systems). However, supercooled liquids can also ex- 
hibit complicated natures, displaying anisotropic mechan- 
ical responses, anisotropic fluctuations [37], and negative 
shear modulus and negative stress correlations. An under- 
standing of the mechanism underlying these complicated 
phenomena will lead us to a better comprehension of the 
violations of the FDT [BEj that they involve, as well as 
certain consequent modifications of the concept of an effec- 
tive temperature [13H21] . Supercooled liquids composed of 
spherical or low-molecular-weight molecules may be excel- 
lent materials for investigating non-equilibrium statistical 
mechanics. 
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Fig. 2. The shear viscosity 77 versus the shear rate y S s at 
various temperatures. The lines are fitted to the function r\ = 
Vso/0- + MTss) + Vf- The arrow indicates the value r/ s o + 77/, 
i.e., the value in the equilibrium state 7 sa = 0. In this study, 
we primarily focused on the sheared non-equilibrium states at 
T — 0.306 and j ss = 10~ 4 (weakly sheared state) and T = 
0.306 and 7 SS = 10 -2 (strongly sheared state), as indicated by 
the black circles. We also considered the different temperature 
case of the state at T — 0.267 and 7 as = 10 -3 , as indicated by 
the black square. 

Table 1. Simulation cases. We performed MD simulations, 
varying the temperature T, the shear rate y S s of the driving 
shear flow, and the amplitude 70 of the oscillating shear strain. 
The item "Figure" indicates the figure that shows the result of 
each case. 
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Fig. 3. The shear moduli (a) G^ and (b) G"H versus the 
frequency u> at T = 0.306 and 7 S3 = 10~ 4 (weakly sheared 
case), with the small amplitude of 70 = 0.01. We also show 
Ge q and G" q in the equilibrium state 7 SS = 0. The symbols and 
lines represent the results of MD simulations and constitutive 
equation ©, respectively. The labels "xy" , "xz" , and "yz" 
denote the component ij = xy, xz, and yz, respectively, and 
the label "eq" denotes the equilibrium state. We indicate the 
two time scales of the constitutive equation (J6]): r s and 17 in 
the sheared state and t s q and Tf in the equilibrium state. 




Fig. 5. The shear moduli (a) G" J and (b) G" SJ versus the 
frequency ui at T = 0.306 and j ss = 10~ 2 (strongly sheared 
case), with the small amplitude of 70 = 0.01. See also the 
caption of Fig. 




Fig. 6. The shear moduli (a) Gl 1 and (b) G"*^ versus the 
frequency uj at T = 0.306 and 7 SS = 10~ 2 (strongly sheared 
case), with the large amplitude of 70 = 0.1. See also the caption 
of Fig. H 



Fig. 4. The shear moduli (a) Gl 3 and (b) G"*^ versus the 
frequency U) at T = 0.306 and j ss — 10~ 4 (weakly sheared 
case), with the large amplitude of 70 = 0.1. See also the caption 
of Fig. [31 
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Fig. 7. The shear moduli (a) G'^ and (b) G"H versus the 
frequency uj at T — 0.267 and y S s = 10~ 3 (different tempera- 
ture case), with the small amplitude of 70 = 0.05. See also the 
caption of Fig. [3] 
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Fig. 9. The shear stress correlation function G(t) in the equi- 
librium state at various temperatures. At long times, the func- 
tion G(t) can be well fitted by the stretch exponential form 
G a exp (— (t/Va) ) as indicated by dotted lines [33], where the 
value r a is the a-relaxation time, and the value G a is the 
plateau modulus. 




Fig. 10. The value 1/(77 — Vf) ~ l/ffao versus the shear rate 
■jss- The temperature is T = 0.267. The straight line is fitted 
by least-squares fit. 
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Fig. 8. The shear moduli (a) Gl 3 and (b) G"*^ versus the 
frequency uj at T = 0.267 and j ss = 10 -3 (different tempera- 
ture case), with the large amplitude of 70 = 0.2. See also the 
caption of Fig. [3] 
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Fig. 12. The shear stress correlation function G 1 ? BB (t). The 
temperature is T — 0.306. The shear rate is (a) j ss = 10 -4 
(weakly sheared state) and (b) 7 SS = 10 -2 (strongly sheared 
state). We also show G eq (t) in the equilibrium state. Again, 
refer to the caption of Fig. for the descriptions of the labels 
"xy" , "xz" , "yz" , and "eq" and the values t s q, t s , and 77. 
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Fig. 13. The response and correlation functions (a) G'? 3 
and G'* 3 T a /ss and (b) G'^ 3 3 and G C H A/ss versus the frequency u 
in the weakly sheared state T = 0.306 and j ss — 10 -4 . We 
also show the functions in the equilibrium state. The symbols 
and lines represent the response functions and the correlation 
functions, respectively. The labels "res" and "cor" correspond 
to the response function and the correlation function, respec- 
tively. Again, refer to the caption of Fig. O for the descriptions 
of the labels "xy" , "xz" , "yz" , and "eq" and the values t s q, t s , 
and Tf. 
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Fig. 15. The effective temperatures (a) T' eS and (b) T" s ver- 
sus the frequency cj in the weakly sheared state T = 0.306 
and 7as = 10~ 4 . We also show the values in the equilibrium 
state, which coincide exactly with T = 0.306. Also refer to the 
caption of Fig. for the descriptions of the labels "xy" , "xz" , 
"yz" , and "eq" and the values t s q , t s , and 77 . 




Fig. 16. The effective temperatures (a) T c ' ff and (b) T'J S versus 
the frequency uj in the strongly sheared state T = 0.306 and 
jss = 10 -2 . See also the caption of Fig. 1151 




Fig. 14. The response and correlation functions (a) G'? 3 s and 
Gcor,-y 3S and (b) G" l3 s and G"™ - versus the frequency u> in 
the strongly sheared state T = 0.306 and jss = 10 -2 . See also 
the caption of Fig. 1131 
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Fig. 17. The susceptibility X^ 33 (t) versus the shear stress 
correlation function G 1 ? {€). The temperature is T = 0.306. 
The shear rate is (a) -y ss — 10~ 4 (weakly sheared state) and 
(b) j ss = 10~ 2 (strongly sheared state). We also show Xeq(i) 
versus G eq (t) in the equilibrium state. The labels u xy" , u xz" , 
and "yz" denote ij = xy, xz, and yz, respectively, and the label 
"eq" denotes the equilibrium state. The slopes are indicated in 
the figures. 



